J 


►  f-n  -Xpcr'y.r  ] 

\ig  vo 


4.  TITLE  AND  SU3TITLE 

Measurements  and  Empirical  Model  Comparisons  of  F-region 
Characteristics  and  Auroral  Oval  Boundaries  During  the 
Solstitial  SUNDIAL  Campaign  of  L987  i 

j  «.  AUTHOR(S)  ,  I 

E.P.  Szuszczewicz ,  P.  Wilkinson,  W.  Swider,  S.  Pulinets, 
j M.A.  Abdu,  E.  Roelof,  T . Fuller-Rowell ,  D.S.  Evans,  T.  | 

i Bateman,  P.  Blanchard,  G.  Gustafson,  R.  Hanbaba,  J.  Joselyn 

it  lf,hH^1e?fiT-wglkfel;t6T’u?rkSS^’  «-«uoheniami 

•7.  PERFORMING  ORGANISATION  NA.MS;3)  V'iC  ACCSSS-vEi)  i  : 

| Phillips  Lab/GPS 

29  Randolph  Road  TP  1 

i  Hanscom  AFB,  MA  01731-3010  II  I  I  l  » 


5-  r'J, 

VDI.VG  NUMBERS 

PE 

61102F 

PR 

2311 

•  TA 

G5 

!  WU 

02 

■!  NO  CRGA.hl 
NUMBiR 


9.  SPONSORING  MC. 


DTIC 

ELECTE  " 


PL-TR-93-2210 


S-.5" '  -ICT:-  AUTHOR  AFFILIATIONS  ON  NEXT  PAGE 
Reprinted  from  Ann.  Geophysicae  11,  601-613  (1993) 


•  12a.  DISTRIBUTION  A  VAIL. \jIL;Tv  S-A  'SMS  I 

;  Approved  for  public  release;  Distribution  unlimited 


13.  435  RACr 


V. .c-rum  <.'.c 


Abstract.  We  report  on  the  SUNDIAL  campaign  con¬ 
ducted  during  the  solstitial  period  May  29-June  7, 
1987.  For  generally  quiet  conditions  a  global  network 
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servations  of  the  Intercosmos  1 806  satellite,  were  com¬ 
pared  with  the  International  Reference  Ionosphere 
(IRI).  The  comparisons  included  F-region  peak  char¬ 
acteristics,  iV„  F,  (i.e.,  /0  F:)  and  hm  F,,  as  well  as  top¬ 
side  height  profiles  .V,  (z).  Overall,  the  IRI  specifica¬ 
tion  of  fa  Ft  was  found  to  be  higher  than  the  observa¬ 
tions.  The  difference  in  the  daytime  hemispheres  was 
nearly  twice  that  in  the  mght-time  hemispheres,  with 
no  obvious  bias  in  universal  time  nor  in  phenomeno¬ 
logical  domains  ;e.g„  equatorial  anomaly,  mid-to-high 


latitudes,  etc.).  With  regard  to  hm  F,  and  topside  pro¬ 
files,  a  small  data  sample  shows  relatively  good  agree¬ 
ment  with  the  IRI.  In  the  topside  tV.  (r)  comparisons 
there  was  good  qualitative  agreement  in  shape;  but 
there  were  quantitative  differences  resulting  from  dif¬ 
ferences  in  :V„  F, .  To  improve  IRI  specifications  at 
high  latitudes  comparisons  also  included  Feldstein, 
DMSP,  and  NOAA/TIROS  auroral  oval  models  and 
their  relative  agreements  with  satellite-bome  particle 
measurements  of  oval  boundaries.  The  NOAA / 
TIROS  model  demonstrates  the  largest  range  in  dy¬ 
namic  oval  boundary  responses  to  magnetic  activity. 
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Abstract.  We  report  on  the  SUNDIAL  campaign  con¬ 
ducted  during  the  solstitial  period  May  29 -June  7, 
1987.  For  generally  quiet  conditions  a  global  network 
of  ionosonde  data,  supported  by  topside  sounder  ob¬ 
servations  of  the  Intercosmos  1806  satellite,  were  com¬ 
pared  with  the  International  Reference  Ionosphere 
(IRI).  The  comparisons  included  F-region  peak  char¬ 
acteristics,  Nm  F-,  (i .e.,f0F2)  and  hm  F2 ,  as  well  as  top¬ 
side  height  profiles  (;).  Overall,  the  IRI  specifica¬ 
tion  of  f0  F2  was  found  to  be  higher  than  the  observa¬ 
tions.  The  difference  in  the  daytime  hemispheres  was 
nearly  twice  that  in  the  night-time  hemispheres,  with 
no  obvious  bias  in  universal  time  nor  in  phenomeno¬ 
logical  domains  (e.g.,  equatorial  anomaly,  mid-to-high 
latitudes,  etc.).  With  regard  to  hm  F2  and  topside  pro¬ 
files,  a  small  data  sample  shows  relatively  good  agree¬ 
ment  with  the  IRI.  In  the  topside  Ne(z)  comparisons 
there  was  good  qualitative  agreement  in  shape;  but 
there  were  quantitative  differences  resulting  from  dif¬ 
ferences  in  Nm  F2 .  To  improve  IRI  specifications  at 
high  latitudes  comparisons  also  included  Feldstein, 
DMSP,  and  NOAA/TIROS  auroral  oval  models  and 
their  relative  agreements  with  satellite-borne  particle 
measurements  of  oval  boundaries.  The  NOAA/ 
TIROS  model  demonstrates  the  largest  range  in  dy¬ 
namic  oval  boundary  responses  to  magnetic  activity. 

Correspondence  to:  E.  P.  Szuszczewicz 
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It  was  also  found  to  be  in  the  best  agreement  with 
observations  and  appears  to  offer  highest  promise  for 
a  proper  high-latitude  adjunct  to  the  IRI  and  an  im¬ 
proved  empirical  specification  of  ionospheric  distribu¬ 
tions  in  auroral  zones. 


Introduction 

O  verview 

Since  October  1984  the  SUNDIAL  Science  Team  has 
conducted  a  series  of  coordinated  measurement  and 
modelling  campaigns  to  understand  better  and  more  ac¬ 
curately  predict  the  global-scale  ionosphere  in  its  quies¬ 
cent  and  disturbed  states.  The  efforts,  which  include  as¬ 
pects  of  solar,  interplanetary,  magnetospheric,  and  ther¬ 
mospheric  physics,  focus  on  understanding  the  coupling 
mechanisms  as  they  manifest  themselves  in  globai-scale 
ionospheric  responses  (Abdu  era/.,  1988  and  1990; 
Leitinger  et  al.,  1988;  Richmond  et  al.,  1988;  Spiro  er  at.. 
1988;  Szuszczewicz  era/.,  1988,  1990,  and  1992;  Biondi 
er  al.  1990;  Emery  et  at..  1990;  Fejer  et  al.  1990;  Kikuchi 
era/.,  1990;  and  Sica  era/.,  1990).  The  effort  has  been 
covering  equinoctial  and  solstitial  periods  separated  by 
9-month  intervals  with  around-the-clock  measurements 
in  8  -  30  day  campaigns  that  have  provided  a  coordinated 
international  data  base  throughout  the  current  solar  cy- 
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cle.  We  report  here  on  the  global  aspects  of  the  third 
campaign,  conducted  during  the  solstitial  period  May 
29-June  7,  1987.  The  results  focus  on  the  testing  of  em¬ 
pirical  model  specifications  of  the  ionosphere  and  the 
auroral  oval  as  the  zero-order  manifestation  of  magneto- 
spheric  inputs;  and  the  measurement/model  results  re¬ 
ported  here  define  the  baseline  conditions  for  other  inves¬ 
tigations  conducted  in  this  campaign  [e.g.  Walker  et  al., 
1991;  Lester  et  al.:  Ruohoniemi  et  al .;  Abdu  et  al.\  Miller 
et  al.:  Denisenko  et  al.  it 993.  all  in  this  issue)]. 


Perspectives  on  program  rationale  and  approach 

The  SUNDIAL  activities  emphasize  the  comparison  of 
observations  with  empirical  and  first-principle  models 
(Szuszczewicz  et  al..  1988;  Schunk  and  Szuszczewicz. 
1988;  Spiro  et  al..  1988;  Sica  el  al.,  1990;  Wilkinson  et  al., 
1992),  in  order  to  provide  further  tests  and  validation  of 
the  models,  improve  their  accuracy  as  appropriate,  and 
advance  our  overall  understanding  of  phenomenologies 
and  associated  cause-effect  relationships.  From  an  em¬ 
pirical  model  perspective,  the  International  Reference 
Ionosphere  (IRI)  has  been  the  model  of  choice  because  it 
is  the  most  widely  used  and  tested  global-scale  ionospher¬ 
ic  model  available  to  date  (Rawer.  1981;  Rawer  and  Ra- 
manamurty.  1985:  Schunk  and  Szuszczewicz,  1988; 
Wilkinson  et  al..  1988).  Our  approach  to  the  IRI  and  its 
comparison  with  our  observations  is  motivated  by  inter¬ 
est  in  it  as  a  baseline  for  time-dependent  studies  and 
determination  of  its  integrity  as  a  global-scale  specifica¬ 
tion  of  conductivity  distributions. 

Relative  to  time-dependent  studies,  one  must  first  ask 
if  a  global-scale  quiet-time  ionospheric  condition  can  be 
defined,  modelled,  and  tested.  This  is  done  more  easily  on 
a  local  or  regional  basis,  but  it  has  yet  to  be  done  on  a 
global  scale.  One  can  point  to  averaged  conditions  (as  in 
the  monthly  averaged  specifications  provided  by  the  em¬ 
pirical  baseline  in  the  IRI).  But  by  its  very  definition,  the 
IRI  is  not  a  "quiet-time"  ionospheric  model,  and  we  sub¬ 
mit  that  the  scientific  community  has  yet  to  identify  and 
test  a  model  that  specifies  the  quiet-time  global-scale 
ionospheric  electron  distribution.  This  is  the  point  from 
which  studies  of  dynamics  should  embark.  And  we  have 
yet  to  answer  the  question  as  to  whether  or  not  the  iono¬ 
sphere  is  ever  in  a  quiescent  state,  that  is,  one  in  which 
variations  are  occurring  only  with  diurnal  and  semi-diur¬ 
nal  variabilities  introduced  in  first  order  by  solar  electro¬ 
magnetic  radiation. 

If  one  can  properly  specify  the  quiet-time  ionosphere, 
then  one  is  in  the  position  to  specify  and  delineate  global 
dynamics  in  ionospheric-thermospheric  interactions  as 
manifested  in  time-dependent  chemical,  electrical,  and 
momentum  coupling  processes  ...  all  of  which  are  depen¬ 
dent  on  the  accurate  specification  of  electron  density  dis¬ 
tributions,  Studies  of  momentum  coupling  and  dynamo- 
driven  fields  involve  the  electron  density  term  in  the  ex¬ 
pression  for  ionospheric  conductivity,  and  there  are  an 
ever-increasing  number  of  large-scale  numerical  models 
that  look  to  empirical  models  like  the  IRI  for  the  specifi¬ 
cation  of  global  conductivity  distributions.  If  that  specifi¬ 


cation  is  biased,  either  on  a  regional  or  global  basis,  or  if 
it  has  diurnal  inaccuracies,  the  large-scale  numerical 
models  will  unknowingly  introduce  these  biases  with  an 
associated  distortion  in  the  results. 

It  is  with  these  perspectives  that  we  approach  the  IRI. 
We  have  developed  a  global  approach  that  integrates  sta¬ 
tion-by-station  results  into  regional  and  global  represen¬ 
tations.  In  this  way  we  look  to  uncover  trends,  biases  and 
inconsistencies,  and  accordingly  point  to  areas  where  cor¬ 
rections  might  be  in  order.  The  work  reported  here  still 
addresses  averaged  representations,  but  in  a  framework 
that  asks  to  what  shorter-period  temporal  baseline  can 
the  monthly-averaged  IRI  be  reasonably  applied.  Here 
we  make  the  comparison  with  a  1 0-day  measurement 
period  and  look  to  biases  in  global  perspectives,  universal 
time,  and  diurnal  characterization.  This  is  done  with  a 
focus  on  the  critical  frequency  (i.e.,  the  density  of  the 
F-region  peak)  as  measured  on  the  bottom-side  by  a 
global  distribution  of  ionosondes. 

Our  bottom-side  model-measurement  comparisons 
are  complemented  by  model-measurement  comparisons 
using  the  Intercosmos  1806  topside  sounder  in  regions 
less  accessible  to  its  ground-based  counterparts.  This  is 
done  as  an  initial  effort  to  test  IRI  specifications  of  top¬ 
side  profiles,  to  expand  the  database,  and  to  establish  a 
cross  check  on  any  trends  that  might  be  observed  in  the 
bottomside  network. 

Finally,  we  address  one  of  the  most  fundamental  defi¬ 
ciencies  in  the  IRI.  its  representation  of  the  high-latitude 
ionosphere,  particularly  as  it  is  affected  by  particle  pre¬ 
cipitation  in  the  auroral  oval  and  by  variabilities  in  the 
position  of  the  oval  in  periods  of  geomagnetic  storms. 
Three  empirical  models  of  the  auroral  oval  (the  Feldstein, 
DMSP  and  NGAA/TIROS  models)  are  reviewed  and 
tested  against  a  data  base  developed  during  the  SUNDI¬ 
AL  period.  This  comparison  establishes  a  quantitative 
comparison  of  the  individual  model  specifications  of  au¬ 
roral  oval  boundaries  (as  a  function  of  geomagnetic  ac¬ 
tivity)  and  provides  the  initial  stage  for  improving  the  IRI 
specification  of  the  high-latitude  ionosphere. 

Our  overall  results  are  cast  in  global-scale  perspectives 
that  include  not  only  the  results  developed  during  the 
1987  solstitial  campaign  but  involve  comparisons  with 
the  equinoctial  study  of  1986.  The  objective  is  the  identi¬ 
fication  of  trends,  biases,  and  regions  for  improvement, 
in  a  way  that  represents  stepwise  progress  toward  an 
accurate  specification  of  the  global-scale  ionosphere  in  its 
quiescent  and  dynamic  states.  The  1986  and  1987  results 
represent  the  solar-minimum  ionospheric  baseline  for  the 
remainder  of  the  SUNDIAL  campaigns  conducted 
throughout  the  ascending  and  solar-maximum  phases  of 
the  current  solar  cycle. 


Global  F-region  measurements  and  mode)  comparisons 

Prevailing  solar-terrestrial  conditions 
during  the  198 /  campaign 

Unlike  the  first  two  SUNDIAL  campaigns,  the  May 
29- June  7.  19S7  campaign  was  solstitial  (as  opposed  to 
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SUNDIAL  CAMPAIGN . 1-10  JUNE  1987 

DATA/MODEL  COMPARISONS:  ALL  STATIONS  ALL  TIMES 


DAY 


PERCENT  DIFFERENCE 


NIGHT 


PERCENT  DIFFERENCE 


Fig.  1.  Distributions  showing  the  number  of  ionosonde  stations 
(the  ordinate)  with  measurements  that  were  larger  or  smaller  (ex¬ 
pressed  as  percent  difference  along  the  abscissa)  than  those  specified 
by  the  1R1  for  the  campaign  period  ofiune  1987.  The  results  were 
integrated  over  ali  times  and  separated  only  according  to  the  day¬ 


time  or  night-time  hemispheres  Percent  differences  at  0.  7.5,  1 5,  30, 
and  60  correspond  to  the  data/ model  comparison  bins  used  in 
Figs.  3  and  4  of  ±5%,  5-10%.  10-20%.  20-40%.  and  >40%, 
respectively.  The  total  number  of  stations  included  in  each  his¬ 
togram  distribution  is  specified  in  the  panels 


equinoctial).  It  occurred  during  the  early  months  of  solar 
cycle  22,  involved  no  major  storm  activity,  and  for  most 
practical  aspects  was  geomagnetically  quiet.  There  were 
no  significant  solar  flares,  the  Covington  index  for  the 
10.7  cm  flux  was  low  (approximately  75  solar  flux  units), 
and  the  solar  wind  characteristics  were  dominated  by 
steady,  low-speed  streams  with  velocities  typically  be¬ 
tween  300  and  400  kms"1,  and  an  excursion  to  speeds 
near  500  km  s  “  1  on  June  2  nd.  The  geomagnetic  field  was 
also  quiet,  with  some  minor  activity  on  May  29  and 
June  6  when  the  3-h  planetary  Kp  indices  approached 
values  near  5.  Otherwise.  Kp  levels  varied  between  0+  and 
2+  throughout  the  other  campaign  days. 


Bottomside  sounder  comparisons 

We  focus  here  on  the  IRI  and  on  ionosonde  observations 
of  peak  F2-region  electron  densities  at  0,  6, 12,  and  18  UT 
covering  the  June  1-10  period  of  1987.  (The  ionosonde 
data  base  was  extended  beyond  the  planned  May  29  to 
June  7  period  to  provide  an  improved  justification  for 
comparison  with  the  monthly-averaged  ionospheric 
specification  of  the  IRI  defined  for  that  month.)  The  re¬ 
sults  of  nearly  50  stations  are  presented  in  two  formats. 
The  first,  depicted  in  Figs.  1  and  2,  integrates  the  model- 
measurement  comparisons  to  look  for  global  trends  iden¬ 
tifiable  with  diurnal  (day/night)  and  universal  time 
(UT  =  0,  6, 1 2,  and  1 8)  variabilities.  The  second,  presented 
in  Figs.  3  and  4,  provides  a  world-wide  view  of  individual 
station  results  with  a  perspective  on  potential  regional 
and  morphological  influences. 

In  the  format  of  a  probability  distribution  of  station 
observations  as  a  function  of  percent  differences  between 
the  measurements  and  the  model,  Fig.  1  integrates  all  the 
resu.s  without  regard  to  region  or  universal  time,  but 
focuses  on  potential  differences  in  day/night  results.  That 
figure  shows  that  the  observation  <uc  skewed  in  the  di¬ 
rection  of  values  less  than  those  of  the  IRI.  Specifically, 


the  comparison  shows  that  the  data  differed  from  the  IRI 
in  the  daytime  hemisphere  by  an  amount  nearly  twice  the 
difference  in  the  night-time  hemisphere  (-7.13%  vs 
—  4.1%,  respectively).  On  an  overall  average,  therefore, 
the  IRI  overpredicted  the  values  of  /0  F,  by  5.6%  (corre¬ 
sponding  to  an  average  overestimate  of  11.2%  in  iVe 
at  the  F-region  peak).  While  these  averaged  values  are  not 
of  alarming  proportions,  the  data  show  a  substantial 
number  of  station  with  the  IRI  specifications  for  /„  F, 
being  15-30%  higher  than  the  actual  measurements.  In 
terms  of  Nt,  this  corresponds  to  a  32-69%  overspecifica- 
tion  by  the  IRI. 

Figure  2  partitions  the  day/night  comparisons  into  the 
UT  =  0,  6,  12,  and  18  periods  to  look  for  trends  in  the 
model-measurement  comparison  with  biases  in  universal 
time.  That  figure  maintains  the  general  trend  of  a  negative 
bias  in  the  observations  when  compared  with  the  IRI, 
without  any  obvious  dependence  on  universal  time.  No¬ 
ticeable  departures  from  this  trend  might  be  attributed  to 
the  UT=0  (daytime)  and  UT  =  6  (night-time)  cases.  The 
former  case  maintains  the  general  trend  (i.e.,  a  negative 
bias  in  the  data  distribution),  but  is  strongly  skewed  with 
the  most  probable  value  being  —15%.  By  comparison, 
the  UT  =  6  (night-time)  case  shows  a  uniform  distribu¬ 
tion,  but  here  the  station  statistics  were  small. 

We  investigate  the  UT  =  0  and  6  cases  more  fully  and 
explore  the  possibility  for  regional  biases  in  the  global 
presentations  in  Figs.  3  and  4  where  the  individual  sta¬ 
tion  results  are  plotted  against  a  background  of  “predict¬ 
ed”  f0  F2  [=8.9(103)  y/N mF2  [cm  “  3J]  contours  from  the 
IRI  for  the  month  of  June  under  near-solar-minimum 
conditions  (SSN  =  30).  The  symbols  (triangles,  circles,  and 
squares)  locate  the  stations  and  quantify  the  level  of 
agreement/disagreement  between  the  observations  and 
the  model  specifications.  (The  quantitative  significance  of 
the  symbols  is  fully  defined  in  the  caption  of  Fig.  3.)  In¬ 
cluded  in  each  plot  is  a  Q  =  4  Feldstein  oval  (Feldstein, 
1963:  Holzworth  and  Meug,  1978)  aUu  vertical  b<us  ai  the 
6  am  and  6  pm  LT  meridians.  (The  Q  =  4  condition  is 
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Fig.  2.  Same  as  Fig.  1  except  the  data  is 
further  partitioned  to  reflect  the  distribu¬ 
tions  at  UT  =  0.  6,  12  and  18  h 


geomagnetically  equivalent  to  Kp  =  2,  an  approximate 
average  for  the  period  under  study.) 

The  daytime  hemisphere  data  at  UT  =  0  (Fig.  3)  is 
dominated  by  stations  in  the  Japan-Australian  meridian, 
stretching  from  sites  inside  the  southern  auroral  oval  to 
northern  latitudes  in  excess  of  60°  geographic.  The  tem¬ 
poral  period  is  primarily  mid-morning,  with  high-,  mid- 
and  equatorial  coverage.  The  data-model  presentation 
shows  that  the  differences  are  everywhere  consistent.  Vir¬ 
tually  without  exception  in  this  meridional  mid-morning 
zone  the  IRI  predicts  values  of  f0  F2  that  are  greater  than 
the  observations,  and  this  is  without  bias  in  morphologi¬ 


cal  domain.  These  same  conclusions  apply  during  after¬ 
noon  hours,  as  can  be  seen  in  Fig.  4. 

Further  inspection  of  Fig.  4  shows  that  station  statis¬ 
tics  are  very  good  (poor)  in  the  daytime  (night-time)  hemi¬ 
sphere.  In  the  daytime  hemisphere,  where  the  trend  has 
been  more  negatively  biased  than  in  its  night-time  coun¬ 
terpart  (for  all  UT),  the  only  departures  from  that  trend  is 
at  high  northern  latitudes  within  Europe  and  the  former 
Soviet  Union.  It  is  these  stations  that  drive  the  most  prob¬ 
able  value  in  the  UT = 6  (daytime)  distnbution  to  0%. 

The  dominant  feature  in  Figs.  1-4  (and  in  UT=12 
and  18  comparisons  of  individual  station  results  not  pre- 
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Fig.  3.  Global  map  of  IRI  contours  (for 
June  with  SSN  =  R  =  30.  UT  =  0)  of  fo  Fz 
(green  contours  in  MHz)  with  Feldstein 
Q  =  4  oval.  Symbols  identify  SUNDIAL 
stations  with  10-day  averaged  faF:  observa¬ 
tions  at  UT  =  0.  J  s  specify  agreement  with¬ 
in  +  5%.  while  □,  (x|.  ,  and  ■  indicate 

that  the  observations  are  less  than  IRI  speci¬ 
fications  by  5-10%,  10-20%,  20  -  40%.  and 
>40%.  respectively.  Q.  g>.  0.  and  #  have 
the  same  quantitative  scaling  but  for  obser¬ 
vations  greater  than  the  IRI  specifications. 
Vertical  lines  identify  the  sunrise  and  or 
sunset  meridian 
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Fig.  4.  Same  as  Fig.  3  but  UT  =  6 


seated  here)  is  that  the  observations  tend  to  be  predomi¬ 
nately  less  than  those  specified  by  the  IRI  without  regard 
to  latitude  or  phenomenological  domain.  This  was  also 
the  case  in  the  SUNDIAL  solar-minimum  equinoctial 
campaign  of  1986.  There  is.  however,  one  point  that  tends 
to  separate  itself  from  the  rest  of  the  1987  data  and  is 
consistent  with  the  1986  observations.  In  both  the  1986 
and  1987  campaigns,  the  observations  over  northern  Eu¬ 
rope  were  consistently  in  better  agreement  with  the  IRI 
than  at  any  other  place  throughout  the  world.  The  con¬ 
clusion  here,  as  with  the  1986  campaign,  points  to  the 
original  formulation  of  the  IRI  which  had  its  largest  iono- 


sonde  data  base  in  the  European  sector.  It  therefore  fol¬ 
lows  that  the  best  agreement  in  fact  should  be  in  that 
domain.  That  continues  to  be  the  SUNDIAL  observa¬ 
tion. 


The  approach  of  averages 

What  appears  above  to  be  a  generally  good  (albeit 
skewed)  agreement  between  model  and  measurements 
should  be  taken  with  some  caution  since  10-day-averaged 
values,  developed  during  four  separate  times  over  a  diur- 
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nal  cycle,  could  provide  unjustified  perspectives  on  the 
model  capabilities  as  a  result  of  undersampled  testing.  We 
provide  some  insight  into  this  perspective  using  full  diur¬ 
nal  characteristics  of  JQ  F,  values  for  Cachoeira  Paulista 
(31  E  long,  —  22.7 '  lat,  —  1 2.4°  mlat)  in  Brazil  and  for 
Canberra  (149  E  long,  —  35. 33  lat,  —43.7°  mlat)  in  Aus¬ 
tralia.  These  data  (in  Figs.  5  and  6,  respectively)  have 
been  selected  since  they  represent  mid-to-equatorial  sites 
separated  by  almost  180°  in  longitude,  and  one  of  them 
(Canberra)  is  included  in  the  Intercosmos  1806  topside 
comparisons  presented  in  the  next  section. 

Figures  3  and  6  show  the  values  specified  by  the  IRI 
for  the  month  of  June  for  a  sunspot  number  equal  to  30. 
The  figures  also  include  two  data  plots.  One  reflects  the 
value  of  f0  F,  averaged  over  the  first  10  days  of  June  (as 
done  in  the  development  of  the  global  comparisons  in 
Fig.  3  and  4)  and  the  second  reflects  the  averaged  values 
for  data  collected  throughout  the  entire  month  of  June. 

The  10-day  averages  are  similar  to  the  30-day  aver¬ 
ages.  This  is  certainly  the  case  in  a  comparison  of  the 
qualitative  features  of  the  diurnal  behaviors  in  both  data 
sets  and  between  the  data  and  the  specifications  of  the 
IRI.  On  a  quantitative  basis,  Canberra’s  10-day  data 
tends  to  be  about  15%  lower  than  the  30-day  data  during 
the  night-time  period  but  in  near  perfect  agreement  in  the 
daytime.  Had  the  30-day  data  been  the  baseline  for  the 
comparisons  in  Figs.  1  and  2,  the  results  would  have  been 
the  same  at  Canberra  for  daytime  measurements  but 
shifted  by  a  5%  -  10%  difference  in  the  positive  direction 
at  night. 

The  plots  for  Cachoeira  Paulista  show  very  good 
agreement  between  the  10-  and  30-day  data  sets,  with 
some  differences  in  the  period  between  12  and  20  UT. 
Both  data  sets  are  in  very  good  agreement  with  the  pre¬ 
sunrise  minimum  but  the  values  differ  markedly  with  the 
sunset  minimum  of  the  IRI.  In  all  cases,  the  10-  and  30- 
day  data  sets  show  values  less  than  those  of  the  IRI.  The 
data  sets  and  the  IRI  also  reflect  the  influences  of  the 
equatorial  anomaly  in  the  late  afternoon  (UT  %  19  hr). 
The  most  significant  difference  between  the  model  and  the 
observations  involves  the  absolute  values  of  f0  Fz,  within 
the  0  and  16  UT  time  frames.  These  differences  are  ap¬ 
proximately  40%  and  30%,  respectively. 


Topside  profiles 

The  accuracy  of  the  IRI  was  further  investigated  by  in¬ 
cluding  comparisons  of  IRI  topside  ionospheric  density 
profiles  with  the  sounder  observations  of  the  Intercosmos 
1806  satellite.  Those  results  are  shown  in  Fig.  7,  together 
with  F-region  heights  hm  Fz  and  peak  densities  Nm  Fz  de¬ 
termined  by  ground-based  ionosondes  at  several  sites  co¬ 
registered  with  the  satellite  soundings.  Figure  7  corre¬ 
sponds  to  data  on  June  4  and  5  in  the  Australian-Japan 
sector  near  150°  E  longitude,  at  an  average  time  of 
UT  =  4.30  h  (LT  =  1 3.8  h): 

Relative  to  a  comparison  between  the  topside  sounder 
and  IRI  profiles  we  offer  the  following  comments:  1.  the 
satellite  data  represent  mid-to-low  latitude  domains  near 
local  noon;  2.  qualitatively,  the  shapes  of  the  Ne  (z)  pro¬ 


files  are  similar,  3.  the  observed  sounder  densities  are 
generally  less  than  those  in  the  IRI  (in  agreement  with 
previously  discussed  results);  and  4.  the  agreement  in  the 
topside  profiles  would  be  improved  if  the  peak  densities  in 
the  IRI  were  reduced  (as  suggested  by  both  the  topside 
and  bottomside  sounder  findings). 

Ground  ionosonde  stations  near  subsatellite  positions 
are  identified  in  the  panels  and  the  plots  include  the  asso¬ 
ciated  ground-station  measurements  of  the  monthly  me¬ 
dian  and  upper  lower  decile  values  for  .Vm  F,  and  hm  F, 
(depicted  as  +)  as  well  as  their  ( Nm  F,,  h„  F,)  instanta¬ 
neous  point  values  at  UT  =  4  h  on  June  4  (depicted  as  •). 
The  topside  and  bottomside  sounder  results  are  seen  to  be 
in  good  agreement.  There  are  several  points  worth  noting 
about  the  results  in  Fig.  7: 

First,  as  a  monthly  averaged  specification  the  IRI  is 
more  properly  compared  with  the  monthly  averaged  bot¬ 
tomside  sounder  determinations  of  Nn  F,  and  h„  F,  (the 
“crosshairs"  in  panels  A-E).  For  all  cases  in  Fig.  7  with 
co-registered  ground-station  data,  the  peak  electron  den¬ 
sities  agree  well  with  the  IRI  and  are  always  within  the 
decile  range  of  the  observations.  The  tendency,  however, 
is  for  the  measured  monthly  median  values  of  Nm  F,  to  be 
less  than  those  specified  by  the  IRI,  consistent  with  the 
results  in  Figs.  1  and  2.  The  peak  height  values  (monthly 
median)  are  also  systematically  lower  than  the  IRI  (apart 
from  Hobart)  but  agreement  is  achieved  within  the  decile. 

Second,  the  IRI  has  often  been  criticized  for  its  topside 
profile  specification  and  its  associated  representation  of 
the  total  electron  content  (see,  e.g.  Bilitza,  1985).  Our 
small  data  sample  suggests  that  the  problem  with  the  IRI 
is  not  so  much  the  shape  of  the  topside  profile  but  rather 
its  reference  point  (,Vm  F, ,  hm  Fz).  With  this  being  a  vari¬ 
able  that  is  locally  adjustable  in  the  IRI.  the  model  agree¬ 
ment  with  observed  topside  profiles  can  be  very  good. 


Empirical  auroral  oval  models  and  observations 

Potentially  more  serious  deficiencies  in  the  IRI  exist  than 
the  inaccuracies  in  tWm  Fz  and  hm  F,.  These  deficiencies 
have  been  discussed  in  earlier  works  (Szuszczewicz  et  al. , 
1988;  Schunk  and  Szuszczewiz.  1988)  where  it  has  been 
pointed  out  that  areas  requiring  attention  include  high 
latitudes  (with  emphasis  on  domains  involving  energetic 
particle  precipitation).  We  have  given  some  attention  to 
the  high  latitude  deficiency  and  executed  a  comparison  of 
empirical  oval  models  that  offer  the  possibility  of  a  mean¬ 
ingful  adjunct  to  the  IRI.  Here  we  discuss  and  compare: 
i.  the  Feldstein  oval  based  on  photographic  data;  2.  the 
DMSP  model  based  on  “in  situ”  particle  data  collected  in 
the  Defense  Meteorological  Satellite  Program  (DMSP); 
and  3.  the  NOAA/TIROS  model  based  on  particle  flux 
and  energy  measurements  conducted  in  a  series  of  NOAA 
environmental  satellite  missions. 


Feldstein  oval 

The  Feldstein  oval  model  (Feldstein,  1963)  is  based  on 
photographic  data  of  auroral  oval  morphology  in  the 


E.  P.  Szuszczevricz  el  ai:  Measurements  and  empirical  model  comparisons  of  F-region  characteristics  and  auroral  oval  boundaries  607 


CANBERRA 


t 

o 

F 

2 

M 

H 


UT 


Fig.  5.  Diurnal  variations  of  /0  F:  at  Canberra,  Australia  ( —  35.32° 
lat.,  149°E.  Long..  —43.63°  dipole  lat.).  The  red  curve  represents  the 
"predictions'-  of  the  IRI.  while  the  green  and  yellow  curves  are  the 
median  ionosonde  observations  for  June  1-10  and  June  1-30. 
respectively 


C  PAULISTA  ( R  - 30) 


Fig.  6.  Diurnal  variations  of  f0  F,  at  Cachoeira  Paulista.  Brazil 
( -22.70;lat..  46.0°E.  long..  -12.40°  dipole  lat.).  The  red  curve 
represents  the  "predictions"  of  the  IRI.  while  the  green  and  yellow 
curves  are  the  median  ionosonde  observations  for  June  1-10  and 
June  I  -30.  respectively 


northern  and  southern  hemispheres  collected  during  the 
International  Geophysical  Year  (from  November  1957 
through  February  1958  in  the  northern  hemisphere  and 
from  May  1958  through  August  1958  in  the  southern 
hemisphere).  The  oval  boundaries  are  represented  by  a 
Fourier  series  (spherical  harmonics)  of  order  4,  with  the 
constants  in  the  series  determined  by  a  least-squares  anal¬ 
ysis  to  minimize  errors  between  data  and  the  fitting  func¬ 
tion  (Holzworth  and  Meng,  1975). 


DMSP  oval 

Data  used  to  develop  the  DMSP  model  (Hardy  et  al. , 
1987)  were  acquired  by  three  satellites,  F2  and  F4  from 
the  Defense  Meteorological  Satellite  Program  (DMSP), 
and  P78-1  from  the  Space  Test  Program.  F4  had  a  circu¬ 
lar  scn-synchronous  orbit  with  an  840  km  altitude  and  an 
inclination  of  97.4°  in  the  1000-2200  h  meridian.  F2.  also 
at  840  km  circular,  was  to  be  sun-synchronous  in  the 
dawn-dusk  meridian  but  its  orbit  precessed  toward  the 
0800  -  2000  h  meridian  over  its  2.5-year  lifetime.  The 
P78-1  satellite  was  in  a  600  km  circular  sun-synchronous 
orbit  in  the  noon-midnight  meridian. 

Data  were  available  from  September  1977 -February 
1980  using  F2,  from  April  1979-August  1980  using  F4, 
and  from  February  1979 -January  1980  using  P78-1.  The 
F2  satellite  recorded  data  from  17  000  northern  and 
southern  hemisphere  passes,  while  10  000  passes  were 
recorded  by  the  F4  satellite.  There  were  1800  hemispheric 
passes  provided  by  the  P78-1  satellite.  However,  only 
spectra  with  the  detector  look-direction  pointed  toward 


the  local  zenith  were  used.  (The  look-direction  of  the  F2 
and  F4  detectors  always  pointed  toward  the  local  zenith.) 

One  complete  energy  spectrum  (50  eV  to  20  keV)  was 
produced  every  second  from  the  F2  and  F4  detectors, 
while  four  energy  spectra  were  produced  every  second 
from  the  P78-1  detector.  To  provide  an  even  distribution 
of  data  over  the  seasons  of  the  year  and  to  provide  suffi¬ 
cient  coverage  at  high  activity,  15  months  were  chosen 
from  the  F2  and  F4  data  and  1 2  months  of  data  were  used 
from  the  P78-1  satellite. 

Using  all  the  spectra  that  fell  within  a  given  MLT  and 
CGMLAT  zone,  the  average  and  standard  deviations  of 
the  differential  number  flux  for  each  of  the  16  detector 
energy  channels  were  calculated.  The  final  result  was  an 
average  spectrum  for  each  zone  at  a  given  level  of  activity. 
These  average  spectra  were  fitted  using  an  Epstein  transi¬ 
tion  function  to  provide  a  representation  of  the  asymmet¬ 
ric  shape  of  the  auroral  oval.  Because  the  coefficients  used 
in  the  mathematical  fitting  procedure  were  slowly-vary¬ 
ing  functions  of  magnetic  iocai  time,  they  were  deter¬ 
mined  by  expanding  their  MLT  dependence  in  a  Fourier 
series  of  order  6  (thereby  reducing  the  number  of  coeffi¬ 
cients  needed). 

In  performing  a  least-squares  analysis  of  the  fitting 
function  to  the  data,  only  values  above  a  specified  mini¬ 
mum  were  used.  This  was  to  exclude  the  values  below  the 
auroral  zone  and  in  the  polar  cap  where  the  four  quanti¬ 
ties  (i.e.,  integral  energy  flux,  integral  number  flux.  Hall 
and  Pedersen  conductivities)  are  roughly  constant  and 
small,  or  zero.  For  the  integral  energy  flux  and  integral 
number  flux,  the  minimum  value  was  106  keV  (cm*  s  sr)~  1 
and  10ft  el  (cirr  s  sr)~ ',  respectively,  to  the  equatorward 


Attitude:  km .  .  Altitude:  km  Altitude:  km 
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Tops'de  Sounder  and  SUNDIAL/1RI  Comparisons 


R»  30.0  Lat.-  -42  00  Long..  144.40 
UT.  June  4,  1987  03:48:24  [13.43  LT] 
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Fig.  7A-G.  Intercosmos  1806  topside 
sounder  profiles  near  UT  =  4.30  h 
(LT  =  13.8  h)  in  the  Australian-Japanese 
sector.  Panels  A-G  correspond  to  sepa¬ 
rate  measurements  {dashed  curves i  with 
the  solid  curve  representing  the  IRI  specifi¬ 
cation.  Observations  A-E  had  subsatellite 
points  in  near-coincidence  with  the  ground 
stations  named  in  the  panels.  The  cross¬ 
hair  symbol  identifies  that  station's  median 
value  for  Nm  F1  and  hmF2,  and  the  upper 
and  lower  deciles  for  the  entire  month  of 
June  at  UT  =  4.  The  solid  dot  identifies  the 
instantaneous  point  value  on  June  4  near 
the  time  of  the  Intercosmos  pass 
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Fig.  8.  Northern  hemisphere  auroral  oval 
mode  comparisons  at  UT  =4.7  for  equiva¬ 
lent  activity  indices.  The  top  panel  shows  the 
Feldstem  oval  with  £)  =  3  The  middle  panel 
shows  the  DMSP  oval  with  K  =  1,  and  the 
bottom  panel  shows  the  NOAA  TIROS  oval 
with  AI  =  3 


side  of  the  maximum  value  and  10'  to  the  poleward  side. 
Likewise,  for  the  height-integrated  Hall  and  Pedersen 
conductivities,  the  minimum  value  was  1  mho  for  both 
the  equatorward  and  poleward  boundaries  (Hardv  et  al.. 
1987). 

The  minimum  height-integrated  Hall  conductivity  of 
1  mho  was  chosen  to  represent  the  oval  boundary  in  our 
comparison  studies  since  Hall  conductivity  involves  both 
number  and  energy  flux,  and  uses  data  from  the  instru¬ 
ment  channels  corresponding  to  energies  higher  than 
about  500  eV  (Robinson  et  al.,  1987).  This  assures  repre¬ 
sentation  of  more  energetic  particle  effects  and  minimizes 
the  contribution  of  "drizzle"  particle  energies  toward  the 
polar  cap. 


NOAA /TIROS  ova/ 

The  NOAA/TIROS  oval  model  (Evans  et  al.,  1988; 
Fuller-Rowell  and  Evans,  1987)  is  based  on  data  from 
the  NOAA/TIROS  weather  satellites  carrying  particle  de¬ 
tectors  capable  of  measuring  energies  in  the  range 
300  eV-20  MeV.  Data  were  obtained  ovc  8-year  peri¬ 
od  from  1978  to  1988  consisting  of  specv.a  from  more 
than  60  000  passes  across  the  polar  regions.  The  satellites 
were  in  sun-synchronous  orbits  at  an  altitude  of  850  km. 

From  a  given  polar  pass,  the  total  precipitating  parti¬ 
cle  power  input  was  estimated  and  assigned  an  activity 
index.  All  passes  within  a  given  activity  index  (AI)  were 
averaged  in  magnetic  latitude  and  magnetic  local  time 
bins.  The  MLT  was  binned  in  8-min  intervals  and  MLAT 
was  divided  into  U  intervals.  Activity  indices  ranged  from 


1  (very  quiet)  to  10  (active).  (More  information  on  the 
binnina  method  and  data  gathering  is  provided  in  Foster 
et  al..  1986). 

The  program  developed  within  the  SUNDIAL  investi¬ 
gation  to  determine  the  NOAA, TIROS  oval  bondaries 
searched  through  the  data  for  entries  with  energy  fluxes 
greater  than  0.25  ergs  (cm2  s)~ '.  The  energy  flux  of  0.25 
ergs  (cm2  s)~  1  was  the  suggested  definition  of  the  aurora 
oval  edge  (D.  E .  ..ns,  private  communication).  By  interpo¬ 
lating  between  the  two  entries,  the  magnetic  latitude  and 
magnetic  local  time  of  the  oval  boundaries  were  deter¬ 
mined  for  a  given  activity  index.  Because  the  boundary 
locations  exhibited  a  “spiky"  structure  when  displayed  on 
a  magnetic  latitude/local  time  plot,  the  oval  edge  data 
were  smoothed  using  a  Fast  Fourier  To  nsform  (FFT). 
This  curve  is  used  to  define  the  auroral  oval. 


Intercomparison  of  the  three  models 

Each  of  the  three  empirical  oval  models  uses  a  differen 
index  to  specify  the  overall  geomagnetic  activity  (and  thus 
the  boundary  of  the  oval).  These  indices  ( Q ,  Kp,  and  AI 
for  the  Feldstem,  DMSP  and  NOAA/TIROS  models,  re¬ 
spectively)  are  all  traceable  to  magnetic  field  fluctuations 
and  are  interrelated  in  a  linear  fashion  as  represented  in 
Table  1. 

To  illustrate  and  compare  the  global  latitudinal  cover¬ 
age  of  these  ova)  prescriptions.  Fig.  8  presents  the  oval 
boundaries  produced  by  each  model  (i.e.,  Feldstein. 
DMSP,  and  NOAA/TIROS,  respectively)  for  a  universal, 
time  (UT)  of  4.7  hr  (corresponding  to  a  midnight  magnet- 
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Table  1.  Relationships  between  the  activity  indices  used  in  the  Feldstein,  DMSP,  and  NOAA  TIROS  auroral  oval  models 
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Fig.  9.  MEAT  of  the  equatorward  boundary  for  a  given  level  of 
magnetic  activity  as  defined  in  the  NOAA/TIROS  ( top  panel),  the 
Feldstein  (middle  panel),  and  the  DMSP  (bottom  panel)  models, 
respectively.  The  bold  solid  (red)  and  dashed  (blue)  lines  identify  the 
extrema  in  magnetic  activity  during  the  May  29-June  7,  1987  peri¬ 
od.  The  triangles  identify  the  upper  and  lower  deciles  of  the  MLAT 
values  in  the  2300x'-75  h  time  domain 


ic  local  time  near  the  O'  magnetic  longitude  meridian)  and 
activity  levels  corresponding  to  Q,  Kp.  and  Al  equal  to  3, 
1,  and  3,  respectively.  These  activity  levels  represent 
equivalent  quiet  conditions  within  the  framework  of  ac¬ 
tivity  indices  defined  in  each  model. 

The  results  in  Fig.  8  are  characteristic  of  our  compari¬ 
sons  of  the  oval  models  for  other  levels  of  magnetic  activ¬ 
ity.  These  characteristics  fall  into  three  general  categories: 
1 )  shape,  2)  oval  thickness,  and  3)  equatorward  extent  for 
comparable  levels  of  magnetic  activity.  Relative  to  shape, 
the  Feldstein  oval  always  represents  a  simple  displaced- 
annula.-ring;  while  DMSP  and  the  NOAA/TIROS  re¬ 
sults  display  a  more  complicated  shape  reflecting  the  na¬ 
ture  of  the  model  fit  to  the  satellite  data.  DMSP  and 
NOAA/TIROS  ovals  are  also  thicker  than  the  Feldstein 
approach  and  they  tend  to  manifest  a  greater  equator- 
ward  penetration  of  the  oval  for  comparable  activity  lev¬ 
els.  The  DMSP  oval  tends  to  reflect  the  most  complicated 
geometry  and  the  most  equatorward  boundary. 

The  dynamic  range  of  the  three  models  and  their  "re¬ 
sponsiveness"  to  magnetic  activity  are  presented  in  Fig.  9. 
The  figure  shows  the  equatorward  boundary  "predicted” 
by  each  of  the  empirical  models  at  2300  ±1.75  h  MLT. 
The  solid  and  dashed  lines  reflect  the  range  in  the  activity 
indices  during  the  SUNDIAL  interval,  alone  with  associ¬ 
ated  positions  of  the  equatorward  boundaries  at  the  ex¬ 
trema  of  observed  activity.  Several  comments  are  in  or¬ 
der:  1)  the  NOAA/TIROS  model  generally  predicts  a 
more  linear  behaviour  and  a  greater  dynamic  range  than 
the  other  two  models;  2)  the  smallest  latitudinal  range 
occurs  in  the  Feldstein  model;  and  3)  the  DMSP  model 
manifests  some  “anomalous”  characteristics  as  a  result  of 
the  mathematical  fitting  function  utilized  in  its  empirical 
formulation. 


Comparison  of  oval  models  with  observations 

To  test  the  relative  accuracy  of  these  oval  prescriptions 
and  illustrate  the  dynamics  of  the  oval  during  the  1987 
campaign  period,  we  focus  on  the  equatorward  boundary 
of  each  of  these  models  near  local  magnetic  midnight  and 
compare  the  results  with  DMSP  satellite  measurements 
of  energetic  particle  fluxes  during  the  May  29  to  June  7, 
1987  period. 

The  satellite  data  were  based  upon  15-s  averages.  w:fh 
a  complete  energy  spectrum  determined  every  second, 
and  the  boundaries  were  identified  as  those  points  at 
which  the  count  rate  for  1-10  keV  electrons  exceeded  the 
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Fig.  10.  Composite  representation  of  Kp  values  and  measured  equa- 
torward  boundaries  of  the  auroral  oval  by  particle  fluxes  on  the 
DMSP  satellite  during  the  period  May  29-June  7,  1987.  The  data 
were  binned  >n  the  time  frame  23  ±  1.75  h  MLT 


Fig.  11.  Comparison  of  the  NOAA/TIROS  model  with  DMSP  par¬ 
ticle  data  observations  of  the  equatorward  boundary  of  the  auroral 
oval  during  the  SUNDIAL-87  campaign 


background  count  rate  by  at  least  a  factor  of  10.  With  this 
definition,  the  measured  equatorward  boundary  of  the 
oval  at  2300+  1.75  h  is  plotted  in  Fig.  10  along  with  the 
prevailing  Kp  values  during  the  SUNDIAL  period.  In¬ 
spection  of  that  figure  shows  that  the  boundary  dynamics 
conform  to  intuition,  moving  more  equatorward  (pole- 
ward)  with  increasing  (decreasing)  values  of  Kp. 

In  comparing  the  data  with  each  of  the  three  models, 
we  found  the  NOAA/TIROS  oval  agreed  best  with  the 
observations  (see  Fig.  11).  Qualitatively,  the  model  and 
the  measurements  show  very  good  correlation  in  tracking 
the  positive  and  negative  excursions  in  the  values  of  the 
Kp  activity  index.  Quantitatively,  the  comparison  reveals 
typical  differences  ranging  from  l°-2°  in  MLAT  and  nev¬ 
er  worse  than  5\  The  biggest  "errors”  occur  during  dy¬ 
namic  transitions  as  on  May  29  and  30  and  again  on 
June  6. 


While  the  comparison  presented  here  is  by  no  means 
exhaustive,  it  reflects  the  merits  of  the  NOAA/TIROS 
oval  model.  Like  the  DMSP  model,  it  is  based  on  "in  situ” 
energetic  particle  data  but  manifests  better  mathematical 
performance  and  includes  an  empirical  specification  for 
the  high  latitude  convection  patterns  that  are  consistent 
with  the  ovals  and  tunable  to  prevailing  geomagnetic 
conditions  (see  e.g.,  Foster  et  al.,  1986).  It  is  also  updat¬ 
able  from  an  on-line  NOAA/TIROS  data  source;  and 
with  the  convection  pattern  capability,  it  lends  itself  to 
adaptive  modeling  techniques  and  a  broader  spectrum  of 
utilities  when  coupled  with  an  ionospheric  model  such  as 
the  IRI. 


Comments,  conclusions  and  future  directions 

Our  results  have  focused  on  continued  testing  of  IRI 
model  specifications  of  the  global-scale  F-region  iono¬ 
sphere  with  complementary  efforts  on  the  potential  for  an 
auroral  oval  adjunct  to  include  the  zero-order  manifesta¬ 
tions  of  magnetospheric  inputs.  Our  approach  to  the  IRI 
and  its  comparison  with  our  observations  of  Nm  F2 ,  hm  F, 
and  topside  profiles  has  been  motivated  by  interest  in  it  as 
a  baseline  for  time-dependent  studies  and  determination 
of  its  integrity  as  a  global-scale  specification  of  conductiv¬ 
ity  distrioutions.  Overall,  we  continue  to  conclude  that 
the  IRI  is  the  best  global-scale  empirical  ionospheric 
model  available  to  date.  We  note  however,  that  the  1987 
solstitial  results  reported  here  as  well  as  the  1 986  equinoc¬ 
tial  SUNDIAL  resuits  showed  that  the  IRI  specifications 
of  /0  F,  were  consistently  higher  than  the  observations. 

We  aiso  note  that  the  empirical  modelling  of  F-region 
characteristics  through  Nm  F,  and  hm  F,  represents  the 
most  simple  component  in  specifying  the  overall  iono¬ 
spheric  morphology.  This  has  dominated  our  approach 
to  date  largely  because  these  parameters  are  generally  a 
routine  data  product  in  ionosonde  observations.  The 
more  difficult  problem  is  faced  at  Fr  and  E-region  do¬ 
mains,  where  the  height  profile  information  can  be  hidden 
by  intense  underlying  layers  and  where  the  density  distri¬ 
butions  can  depart  in  major  ways  from  the  classical  lam¬ 
inar  descriptions  of  the  E-  and  F^domains.  This  is  partic¬ 
ularly  true  in  the  presence  of  intermediate,  transitional, 
descending,  or  otherwise  sporadic  layers  (Wilkinson 
et  al .,  1992).  These  layers  are  produced  by  thermospheric 
winds,  electric  fields,  and  ion  composition  effects,  and 
they  can  play  an  important  role  in  modifications  of  lower 
ionospheric  conductivities  and  dynamo-driven  fields.  It  is 
in  this  lower  ionospheric-thermospheric  domain  that  ma¬ 
jor  advances  in  empirical  modeling  are  important.  It  is 
therefore  necessary  that  the  diurnal,  seasonal,  and  solar- 
cycle  variability  of  these  intermediate  and  sporadic  layers 
be  properly  defined,  not  only  to  upgrade  empirical  model 
specifications  but  to  provide  an  empirical  baseline  upon 
which  to  test,  validate,  and  improve  first-principle  mod¬ 
els.  This  type  of  effort  has  been  started  recently  in  the 
SUNDIAL  program,  with  some  advances  made  in  the 
diagnostics  and  analysis  of  these  layers  through  careful 
treatment  of  ionogram  data  coupled  with  first-principles 
modeling  support  (Wilkinson  et  al.,  1992).  This  effort  will 
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continue  on  a  global  scale,  with  the  intention  being  an 
improved  empirical  specification  of  this  lower  ionospher¬ 
ic  domain. 

This  lower  region  and  its  associated  plasma  structures 
are  also  modified  in  important  ways  by  energetic  particle 
precipitation  at  high  latitudes.  This  has  long  been  a  criti¬ 
cism  of  the  IRI.  Our  initial  efforts  to  mitigate  this  problem 
has  been  to  explore  the  potential  to  include  an  empirical 
specification  of  the  position  and  dynamics  of  the  auroral 
oval  boundaries.  At  a  minimum,  this  will  help  identify  the 
location  and  movement  of  the  oval  (at  least  as  made  pos¬ 
sible  by  the  averaging  processes  built  into  the  models)  and 
facilitate  the  identification  of  cause-effect  terms  for  relat¬ 
ed  aspect  of  global-scale  ionospheric  dynamics.  Future 
efforts  suggest  that  associated  enhancements  in  the  iono¬ 
spheric  densities  should  also  be  incorporated  into  the  IRI 
in  order  to  reflect  ionization  processes  resulting  from  the 
precipitating  energetic  particles.  The  appropriate  strategy 
here  is  not  clear,  due  in  part  to  the  complex  nature  oi  the 
spatial  and  temporal  variability  of  particle  precipitation 
patterns.  Future  activities  should  investigate  this  aspect 
of  the  empirical  modelling  problem,  with  the  intention 
primarily  focused  on  the  use  of  such  obvious  inputs  as  the 
energy  and  density  fluxes  specified  in  the  DMSP  and 
NOAA  TIROS  models.  Referring  now  to  our  compari¬ 
sons  of  these  two  models,  as  well  as  the  much-used 
Feldstein  oval,  we  note  that  our  tendency  is  to  favor  the 
approach  of  the  NOAA/TIROS  results,  partly  because  of 
its  (apparently  proper)  dynamic  range,  its  inclusion  of 
high-latitude  convection  field  patterns  (derived  self-con- 
sistently  with  the  oval  boundaries  and  cross  polar-cap 
potentials),  and  the  model  agreement  with  the  observa¬ 
tions  conducted  during  the  June  1987  campaign. 

Primary  deficiencies  in  the  IRI  continue  to  be  its  repre¬ 
sentation  of  high-latitude  phenomenology  and  the  speci¬ 
fication  of  E-  and  F-region  densities  in  a  way  that  allows 
the  accurate  inclusion  of  the  E/F  region  valleys  as  well  as 
intermediate,  descending,  and  transitional  layers.  Ad¬ 
vances  are  now  being  made  in  those  areas,  and  it  can  be 
expected  that  there  will  be  improvements  in  the  IRI,  its 
representation  of  monthly-averaged  ionospheric  proper¬ 
ties,  and  in  its  adaptability  to  a  time-dependent  descrip¬ 
tion  of  dynamic  events. 
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